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Abstract. We have investigated the problem of the distribution of both masses and orbital radii of planets resulting from the 
gas-accretion, gas-capture model. First we followed the evolution of gas and solids from the moment where all solids are in the 
form of small grains to the stage when most of them are in the form of planetesimals for a set of different initial masses and 
sizes of protoplanetary disks. Based on that we performed Monte Carlo calculations describing the formation of giant places 
at different locations. We included the effects of type II migration and growth of the mass of the planet after the gap opened. 
We discuss how these effects influence the final distribution of giant planets. We show that when the giant planets are not able 
to migrate or grow in mass after the gap opens, their distribution is mainly determined by the properties of the gaseous disk. 
However, with those two effects included, reproducing the parameters of the gaseous disks from the distribution of planets 
becomes difficult. We also checked the roles of both the material of which the solids consist and the mass of the central star. 
The main result is that, in disks around less massive stars, giant planets at the given location tend to be less massive. At the 
same time, the giant planets with the given mass tend to form closer to the less massive stars. 



1. Introduction 

Current radial velocity surveys have led to the discovery of 
over 150 extrasolar planets around m ain-sequence stars (see 
Mar cv et alJl2005t iMavoret al.l Eo04). The large majority of 
these planets are probably gas giant planets, as their masses 
are above 100M e . Such a collection provides a good data set 
for c omparing predictio ns of theories of planet f ormation (e.g. 



Ilda & Linl2004albHAlibert et al.l2005tlkornet et alJ2005h . 

The standard model for the formation of giant planets is 
the core-accretion, gas capture-model. The numerical calcula- 
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tions of Pollack et al. ( 1996) show, that the formation of a giant 
planet in this model can be divided into three phases. During 
the first one, the solid core is formed by the collisional accumu- 
lation of planetesimals. The second phase begins when the core 
reaches the mass of a few Earth masses and starts to accrete a 
significant amount of gas. During this phase the envelope stays 
in quasi - static and thermal equilibrium, as the energy radi- 
ated by the envelope is compensated for by the energy released 
by accreted planetesimals. As during this phase the protoplanet 
accretes gas with a higher rate than solids, the mass of the en- 
velope finally becomes equal to the mass of the core. At this 
moment phase 3 begins, during which the planet rapidly grows 
in mass by runaway accretion of gas. The final mass and loca- 
tion of a giant planet are determined by its gravitational interac- 
tion with its environment. As it grows in mass it induces spiral 
waves in the gaseous disk. This leads to the transfer of angular 
momentum resulting in the inward migration of the planet and 



possibly in the gap opening dLin & Papaloizoull986tlLin et alJ 
1996; W ardll 19971) . This last phenomenon strongly reduces the 
further growth of the planet. 

The main problem with that scenario is related to the 
timescale required to form a giant planet in it. In general it is 
the same order of magnitude as the lifetime of the protoplan- 
etary disk, and it is not necessarily certain if the giant planet 
is able to form before dispersion of the disk. Close to the star, 
the formation time of a giant planet in the gas capture model 
is determined by phase 2, while at larger distances (> 10AU) 
the lengths of phases 1 and 2 become comparable. The lengths 
of these two phases depend on the initial surface density of the 
planetesimal swarm in the given location. The larger the den- 
sity, the faster the core grows and reaches higher mass at the 
end of phase 1 . With the higher mass of the core, the length of 
phase 2 also diminishes. In general, at every distance from the 
star there is a critical value of the surface density of planetesi- 
mals which enables formation of a giant planet before disper- 
sion of the protoplanet ary disk; for a more detailed discussion 
see lKornet et al.l l l2006l) . 

However, the density of the protoplanetary swarm is not 
in a simple relation with the density of gas in the disk from 
which it emerges. While the gaseous component evolves in a 
viscous way, the evolution of the solid component is governed 
by proce sses of gas drag, coagulation, s edimentation, and evap- 
oration dWeidenschilling & Cuzzill993h . A significant redistri- 
bution of solid material takes place in the result, and in the 
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inner parts of the disk its surface density can be significantly 
enhan ced compared to the in itial value (Stepinski & Valageas 
l!997t IWeidenschillir^l2003l) . Consequently, analysis of the 
possible masses and orbits of giant planets resulting from the 
core accretion scenario should also include the global evolution 
of solids in protoplan etary disks. A simpl e model of this evolu- 
tio n was propos ed by Kornet et al. | l|2 0Qll) and f urther extended 
bv lKornet etail J2002I) and llCornet et all l2005h to include sub- 
sequent formation of giant planets. 

In this paper we extend our models to also include the ef- 
fects of migration and the gap opening by the planet. It enables 
us for the first time to characterize the distribution of plane- 
tary masses and final locations resulting from the core accre- 
tion model, not only including the growth of the planet but also 
the preceding evolution of solids, which determines the surface 
density of the planetesimal swarm. In Sect. [2] we explain our 
approach to the evolution of the protoplanetary disk and planet 
formation. The results are presented in Sect.[3]and discussed in 
Sect. E| 

2. Methods 

In our approach we divide the formation of giant planets into 
three phases in a natural way. In the first one the planetesimal 
swarm is formed in the protoplanetary disk by the collisional 
accumulations of solids. This phase lasts till the gravitational 
interactions between solid bodies become the dominant factor 
governing their evolution. In the second epoch the giant planet 
is build by the accretion of the planetesimals, and at later stages 
of gas onto the protoplanetary core. Finally, when the mass of 
the planet reaches a critical value, the planet migrates within 
the disk and reaches its final position. We now discuss how we 
modeled each of these phases. 



2.1. Formation of a planetesimal swarm 

We describe the protoplanetary disk as a two-component 
fluid consisting of gas and solids. We modeled the 
gaseous component as a g eometrically thin turbulent a disk 
(Shakura& S unvaevlll973l) . Its surface density E is given as 
a function of distance a from the star an d time t in terms of 
a selfsimilar solution of Stepinski All other quantities 

characterizing the gas were obtained by solving t he stan dard 
set of equations for a thin-disk approximation (e.g. lFrank et alJ 
The initial conditions were parameterized by two quan- 
tities: the initial mass of the disk Mo and its initial outer radius 
So. 

The crucial approximation underlying our approach to the 
evolution of solids is that the size distribution of particles at 
any given radial location is narrowly peaked around a mean 
value particular to this location and time. In the practical imple- 
mentations it means that sizes of the particles s are expressed 
as a single-value function of time and position, s = s(t, r). 
Collectively, all solid particles are treated as a turbulent fluid 
characterized by its mean surface density E s (f, r), which is gov- 
erned by the continuity equation: 

3E S 1 d 

-^ + -^-(rV s E 8 ) = 0. (1) 
ot r or 



The radial drift velocity of solids V s is the result of the frictional 
coupling between solids and the gas and depends on the local 
size of the particles and properties of the gas. In the regions 
where the temperature of the disk is higher than the evapora- 
tion temperature r evap , the solids are treated as vapour with a 
velocity equal to the velocity of the gas. 

The size of the particles is governed by the second equation: 
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7^ + -|-(>V s E size ) = /E s 
ot r or 



(2) 



where E s j ze = sE s . The source function / describes the growth 
of particles due to mutual collisions and subsequent coagula- 
tion. The main assumptions used in its derivation are that all 
collisions between particles lead to coagulation and that the rel- 
ative velocities of collid ing particles are g i ven by the turbulent 
model described by Stepinski & Valageas ( 1997). In the calcu- 
lation of the density of solids from their surface density, the 
effect of their sedimentation toward the midplane of the disk 
is taken into account by evolving the scale height of the solid 
disk. 

Initially the surface density of solids is everywhere a con- 
stant fraction of the surface density of gas, and their sizes 
amount to 10~ 3 cm. The results discussed in the subsequent 
sections do not depend on the choice of that particular value, as 
long as the solids are initially small enough to couple well to 
the gas. 

2.2. Growth of planets 

Our procedure for the formation of the giant planet starts 
when the planetesimals at a given point in the disk reach 
radii of 2 km. That moment determines the initial surface 
density of the planetesimal swarm. The growth of the proto- 
planetary core is described b y the following formula given by 
Papaloizou & Terauem ( 1999): 



M c - CiC cap R p R}iQ.KT's 



where 



Rh = a 



3M* 



1/3 



(3) 



(4) 



is the radius of the Hill sphere of the planet and Qk 
the K eplerian angular veloci ty. The value of C\ given by 
Papa loizou & TerquerrJ dl999l) is 81^/32; we use a factor of 
5 (the difference comes from the different definition of R^). 
The quantity C cap describes the increase in the effective cap- 
ture radius of the planet within respect to its real radius R v 
due to the interaction of planet esimals with the envelope of 
the planet (tPodolak e^ajjjlgg sfc. We use for it an approxi- 
mate fit provided by IHubickvil d20o"il) . made to the results of 
Bodenheimer et al. (2000). For those core masses that are less 
then 5M e , no increase in the effective capture radius is as- 
sumed, i.e. C cap (M c < 5M e ) = 1. For larger core masses, it 
increases linearly with the mass of the core, reaching its max- 
imum value of C cap = 5 for M c = 15M e . We assume that the 
surface density of planetesimals Ed is constant throughout the 
feeding zone which, extends to 4 Hill radii on both sides of the 
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planetary orbit. Furthermore, changes only due to accretion 
of the planetesimals onto the core and the growth of the feeding 
zone. 

Significant accretion of gas onto the core begins when it 
reaches a critical mass of 



M. 



= 10 



M c 



0.25 
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^env j 




llcm* g 1 j 



0.25 



(5) 



\Q- b M @ yr 

Jlkoma et alJ20 00). where /c env is the opacity in the envelope of 
the planet. Its actual magnitude is currently poorly constrained. 
We assume that /c env = 1 cm 2 g _1 . When the mass of the pro- 
toplanet exceeds M CiCI j t , it contracts on t he Kelvin-He l mholtz 
time scale tkh- By fitting the result of IPollack et alJ \ 199d) 
iBrvden et alT (EoOO) show that 
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TKH^lO^f^, , , , , 

\ M m / \ 1cm 2 g" 1 



yr 



(6) 



where the exponents b — 10 and c — 3. However, the gas accre- 
tion rate connected with this contraction cannot be higher than 
the speed with which viscous transport of gas replenishes the 
feeding zone of the planet. Consequently, we adopt the follow- 
ing equations for the gas accretion rate onto the star 



dM P 



df 



(7) 



where Mii s k is the accretion rate in the disk without the planet. 
We assume that it is equal to naC s lTL, where Cs is the sound 
speed in the gaseous disk, H its scale height, and S its gas sur- 
face density at the given location. 

When the Hill radius of the planet becomes larger than 1 .5 
of the scale height of the disk, the planet induces a strong tidal 
torque on the disk and opens a gap in it. As a result the accretion 
of gas is strongly diminished. Its maximum rate is then equal 
to: 



M env = M disk \ 1 .668 | — | exp 



Mr 



1/3 



l.5M 3 



+ 0.04 



(8) 



where Mj is the Jupiter mass JVeras & Armitag el20041. 



2.3. Migration 

The gravitational interaction of a planet with the disk leads 

to its migration and to formation of the gap in the disk 

I * ' II I I II —i 

( Lin & Papaloizou 1986; Ward 1997). For a low-mass planet 
the interaction is linear and results in so-called type I migra- 
tion. In contrast to this scenario, high-mass planets open a gap 
in the disk that reduces the timescale of migration (referred to 
as type II migration). 

Analytical calculations show that the timescale of type I 
migration in a laminar disk is much shorter than both the disk 
lifetime and the timescale of planet formation. Consequently, 
if type I migration has taken place, nearly all planets would 
be accreted onto the star and the probability of finding any 
planetary system would be very low. Likewise, recent simula- 
tions bv lNelson & Papaloizou! J2004I) suggest that the torques 
exerted onto a low mass planet in a turbulent MHD disk fluctu- 
ate strongly. As a result the planet proceeds at a random walk 



and the mean value of the migration velocity is highly reduced. 
Due to these uncertainties we neglect the type I migration in 
our models. 

Type II migration begins when the planet is massive enough 
to open a gap in the gas. It happens when the Hill radius of the 
planet is larger than 1.5 of the disk scale height. If the mass of 
the planet is negligible in comparison to the mass of the disk, 
the inward velocity of the planet is regulated by the viscosity 
of the disk llWardl 19971) : 

df al\ 

where C s is the sound velocity in the gas. In the case of a planet 
with a larger mass than the mass of the disk, the migration is 
slowed down. The variation in the planetary angular momen- 
tum is then equal to the angular momentum transport rate in 
the disk (Lin et alJl996l) : 

^-[M p fl 2 Q k ] = -^vEQ k fl 2 . (10) 
at 2 



In our calculations we used the smaller of the values of da/dt 
given by Eqs. and JlOi . The migration of the planet was 
calculated using the 4th order Runge - Kutta method. It stops 
when either the time from the beginning of the calculation (in- 
cluding the time needed for the formation of the planetesimal 
swarm) is longer then the lifetime of the protoplanetary disk 
Tf, or the orbit of the planet becomes smaller than 0.01 AU. In 
the latter case we assumed that the planet is accreted onto the 
star. In our calculation s tv — 3 X 10 6 yrs, which agrees with 
observations (Haisch et al. 2001). 

3. Results 

Using the methods described in the previous section we were 
able to investigate the range of masses and orbits of giant plan- 
ets which are able to form in protoplanetary disks character- 
ized by different initial conditions. First, we calculated a grid 
of models of protoplanetary disks. Every model is character- 
ized by two parameters: the initial mass (Mo) and outer radius 
(Rq) of the gaseous disk. Their initial masses are uniformly dis- 
tributed between 0.02 and 0.2 M . The range of outer radii was 
chosen with the goal of including all models in which the for- 
mation of giant planets is possible. The number of radii is con- 
stant per interval of lg.Ro- The viscosity parameter a is equal 
to 0.001 for all models. The solids consist of water ice with a 
bulk density of 1 g cirT 3 and have an evaporation temperature 
of 150 K. For every value of [Mo, Rq] we also checked the grav- 
itational stability of the corresponding gaseous disk. In some 
cases the value of the Toomre parameter 



Q 



(11) 



drops below 1 in the outer regions of the disk. It means that 
they are gravitationally unstable with respect to the axisym- 
metric modes. We assumed that these parts of the disk fragment 
and some giant planets are formed there on a very short time 
scale (see i.e. |B2ssl2002). Consequently we used modified ini- 
tial values of the mass and outer radius of the disk, which cor- 
respond to the mass and size of the stable part of the original 
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disk. Our distribution the of initial parameters of protoplanetary 
disks is an artificial one. In general one should use the distribu- 
tion that occurs in nature. Unfortunately current observations 
do not provide such information. However, our results remain 
valid, as we are interested only in characterizing the general set 
of possible orbital sizes and masses of planets. 

Using that grid of evolved models we performed Monte- 
Carlo calculations to produce the M p - a distribution of planets 
resulting from our models. With the uniform probability we 
randomly chose one of the models from the calculated grid. 
Next, we chose the initial distance of the planet from the star 
a with constant probability per interval of lg a. We evolved the 
mass and semimajor axis of the planet according to the equa- 
tions given in Sects. I2~2l and l2. 31 starting from the time at which 
the planetesimals at radial distance a reach a radius of 2km in 
the given model of the protoplanetary disk. To better illustrate 
how the type II migration and the accretion of gas after opening 
the gap change the parameters of the giant planets we also per- 
formed calculations in which those two phenomena were not 
included. 

The results of calculations in which the planets were nei- 
ther permitted to increase their masses after the gap opening in 
the disk or to change their semimajor axes are presented in the 
upper panel of Fig.^ It shows that the planets can be divided 
into two main groups (see also Fig. [2}. In the first one there are 
bodies with masses smaller than ~ 20M ffl . These bodies were 
not able to accrete a significant amount of gas so the nearly 
only consist of their solid cores. In the second group, there are 
planets with masses larger than ~ 100M e . In these cases the 
runaway gas accretion has already taken place and the majority 
of their masses are in the form of gaseous envelopes. They are 
massive enough to open a gap in the disk. The small number 
of planets with intermediate masses results from the fact that, 
as the mass of the planet grows above ~ 20M e , the timescale 
of the accretion of gas quickly becomes much shorter than the 
lifetime of the protoplanetary disk. Consequently, the probabil- 
ity of the disk to disperse be fore the planet gro ws sufficiently 
to open a gap is very low (see llda & Linl2004al) . 

The giant planets in this set of models form in a region 
between 1 and 20 AU from the central star. They can be divided 
into two subgroups, each of them forming one of the wings 
of the V-like shape of the distribution of giant planets on the 
a - M p diagram. As the final mass of the giant planets in this 
set of models is determined by equation 

(1 5hV^ 
~a~) M * ' (12) 

those two subgroups of planets reflect different ways in which 
the HI a ratio depends on distance from the star a at the moment 
of gap opening. Giant planets at distances larger than ~ 10AU 
are at that time in the opticall y thin parts of th e planetary disk 
in which (H/a) ~ a 1 ^ 20 (see Stepinskil ll998l) . Consequently, 
the masses of these planets are correlated with their distances 
from the star. The correlation is extremely narrow because scale 
height of the disk does not depend on its initial mass and outer 
radius in this regime in the used model of the gaseous disk. 
Giant planets at orbits smaller than ~ 10 AU in the moment of 
gap opening are in the optically thick part of the disk in which 
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Fig. 1. The distribution of theoretically predicted semimajor 
axes and masses of planets around a star with a mass of 1M . 
The upper panel shows planets that neither migrate nor grow 
in mass after they are massive enough to open the gap in the 
disk. The lower one presents results for planets that undergo 
type II migration, but do not grow after gap opening. The solid 
material in the disk is water ice. 

(H/a) ~ aT llA , so in their case their masses are anticorrelated 
with the distance from the star. The spread of this correlation 
is larger, because in that case the exact form of the H(a) rela- 
tion depends on the initial mass and radius of the underlying 
disk. This subgroup contains the most massive, closest plan- 
ets formed in this set of simulations. They have masses around 
450M ffi and orbits around 4 AU. 

In the next set of models we include the effects of type II 
migration, while planets still do not accrete any more after gap 
opening. The dependence of the masses of the planets resulting 
from this simulations on the sizes of their orbits is presented in 
the lower panel of Fig.^ As expected, the only differences in 
comparison with the previous case are in the region of the giant 
planets. Moreover, the range of the masses of these planets in 
both figures are the same. Giant planets at larger orbits open the 
gap later because the process of their formation is longer due 
to lower densities of planetesimal swarm and smaller Keplerian 
frequencies (see Eq.[3}. Consequently the changes in their or- 
bits are smaller as they have less time to migrate. Also the ac- 
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Orbital radius 

Fig. 2. The schematic diagram of distribution of theoretically 
predicted semimajor axes and masses of planets around a star 
with a mass of 1M without effects of the gap opening and 
subsequent accretion of gas onto the planet. 

cretion rate in the disk that determines the speed of migration 
is lower at later times. In our calculations, only planets at or- 
bits smaller than ~ 9 AU are able to shrink their orbits by a 
factor larger than 10%. As the result, the distribution of plan- 
ets on the a - M p diagram for a > 9 AU is nearly the same as 
in the case without migration. In the case of planets formed at 
smaller distances from the star, some of them are able to mi- 
grate to orbits as small as ~ 2AU. Additionally, more massive 
planets tend to shrink their orbits by a larger factor. This can be 
explained by the higher rate of migration in the hotter disks. At 
the same time such disks have larger scale heights, and giant 
planets can reach larger masses before gap opening and stop- 
ping their growth. All these factors lead to the small spread of 
giant planets around a curve on the a - M p diagram. 

In the third set of calculations, the planets both migrate and 
grow in mass after gap opening. The distribution of the result- 
ing planets is shown in Fig. [3] The range of their final semima- 
jor axes is similar to the one in the previous set of models. The 
maximum size of planets at the given distance is ~ 2000M e at 
2 AU and decreases to ~ 400M e at 17 AU. The increase in the 
masses of the planets, after gap opening tends to be anticorre- 
lated with the semimajor axes of their orbits. We explain this 
fact as follows. The rate at which a giant planet accretes gas at 
this stage of its evolution is determined by the accretion rate 
in the disk, and as such does not strongly depend on the dis- 
tance from the star. At the same time the protoplanetary cores 
at smaller distances grow on shorter timescales. There are two 
reasons for this. First, the accretion rate of planetesimals onto 
the core is proportional to Keplerian frequency; second, the 
surface density of the nascent planetesimal swarm tends to be 
larger closer to the star. Consequently, protoplanets at smaller 
orbits are able to start accreting significant amounts of gas for 
a longer time before dispersion of the disk. Additionally, the 
accretion rate in the disk, which determines the accretion rate 
of gas onto the planet, is larger in this earlier epoch. As a re- 
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Fig. 3. Effect of planet growth after gap opening on the theoret- 
ically predicted distribution of orbital sizes and masses of plan- 
ets around a star with mass of 1M . The upper panel shows the 
results for planets that undergo type II migration and growth 
after gap opening. The lower one shows results for planets that 
undergo type II migration, but not grow after gap opening. The 
solid material is water ice. 

suit, the maximum size of the planet is everywhere a decreas- 
ing function of the distance from the star, even in those regions 
where the relation was the opposite at the moment of gap open- 
ing. 

Next, we investigate how the relation between the mass of 
the planet and semimajor axes of its orbit depends on the sort of 
material solids consist of. For this purpose we performed sim- 
ilar simulations to those described above, but with solids made 
of high temperature silicates instead of water ice. For the evap- 
oration temperature of silicates we adopted a value of 1350 K, 
and for their bulk density 3.3 g cm 3 . The upper panel of Fig.|4] 
presents the distribution of planets resulting from these calcu- 
lations. Because the solids can survive at higher temperatures 
than ice, the giant planets in that case can form closer to the star. 
Because the rate of migration does not depend on the composi- 
tion of solids, their final orbits can also be smaller. At the same 
time, the distribution of the final surface densities of the plan- 
etesimal swarms, as a function of the distance from the star, 
seems to be a natural extension of this distribution in the case 
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a [All] 

Fig. 4. Effect of different kinds of material that the solids in 
the protoplanetary disk consist of on the theoretically predicted 
distribution of orbital sizes and masses of planets around a star 
with a mass of 1M Q . The upper panel shows the results from 
models in which solids are composed of silicates, the lower 
one - from water ice. In both cases the planets undergo type II 
migration and grow after gap opening. 



of water ice solids. As a result, the relation of a - M p for planets 
with silicate cores at orbits smaller then 6 AU is just an exten- 
sion to smaller values of a of the same relation for planet with 
water ice cores. 

Consequently, in our models the most massive planets (with 
M p ~ 5000M ffi ) tend to be located closest to the star at or- 
bits with semimajor axes ~ 0.5 AU. On the other hand, plan- 
ets with silicate cores at orbits larger than 6 AU tend to have 
smaller masses than their ice core counterparts. The reason is 
that at these distances the time during which the core of the 
protoplanet only accretes solid material becomes a significant 
fraction of the whole time of its formation. The rate of the ac- 
cretion of planetesimals onto the core with the same mass, but 
composed of more dense material, is lower due to its smaller 
physical dimensions. As a result, planets with silicate cores 
start their gas accretion later and grow to smaller masses be- 
fore the dispersion of the disk. At distances larger than ~ 8 
AU, we do not obtain any giant planets with silicate cores. 




a[AU] 



Fig. 5. The distribution of theoretically predicted semimajor 
axes and masses of planets around stars with masses of 1M 
(lower panel) and O.5M (upper panel). In both cases the plan- 
ets undergo type II migration and grow after gap opening. 

Finally, we investigated the influence of the central star 
mass on the relation between masses and orbit sizes of giant 
planets. In Fig. |5]we present the results for two values of the 
stellar mass: 0.5 and 1 M Q . The planets around less massive 
stars also tend to be less massive and have slightly smaller or- 
bits. The first fact is explained by the lower accretion rates in 
these disks, as these rates also determine the maximum rate at 
which the planet can accumulate gas. The tendency of giant 
planets to form around less massive stars at smaller radii is the 
result of (a) a more effective radial redistribution of solids in 
the disk prior to the formation of planetesimals and, as a re- 
sult the higher surface density of the planetesimal swarm, (b) a 
lower value of the minimal surface density for the formation of 
a giant planet around less massi ve stars at orbits smaller than 
~ 10 AU (see lKornet et all2006l) . 

4. Conclusions 

We have presented the results of simulations describing the for- 
mation of giant planets. Our model include all important phases 
of the core-accretion scenario, beginning from the formation of 
the planetesimal swarm from smaller solids, until the gap open- 
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ing in the disk and subsequent planet migration. We presented 
the distribution of a - M p for planets resulting from our models 
for a set of different initial masses and sizes of the protoplan- 
etary disk and for the place of planet formation. However, our 
results cannot be compared in a simple way to data for extra- 
solar planets, because of the lack of knowledge of the distribu- 
tion of the initial parameters of protoplanetary disks that occur 
in nature. Instead our aim was to characterize the relationship 
between the final orbital radii and masses of planets that are 
in general possible to obtain and to check how the resulting 
a - M p distribution changes when including different physical 
processes such as migration and gas accretion by the planet af- 
ter gap opening. We also investigated whether the distribution 
of giant planets is mainly determined by the parameters of the 
gaseous disk or the distribution of the planetesimal swarm. 

We have shown that, if we do not include the effects of mi- 
gration nor accretion of gas by the planet after gap opening, the 
a—M v distribution for giant planets is mainly determined by the 
parameters of the gaseous disk. In that case one is able to read 
the information about the dependence of the scale height of the 
gaseous disk on the distance from the star from the distribution 
of giant planets. The minimal and maximal distance from the 
star gives information about the surface density of planetesimal 
swarm at those places. 

By including migration in our models, reconstructing the 
the properties of the gaseous disk from the distribution of gi- 
ant planets becomes more difficult. The distance that the planet 
can move inward depends not only on the gaseous environment 
but also on the time between gap opening and dispersion of the 
disk. But the first of these two moments is determined by the 
surface density of planetesimal swarm. This, on the other hand, 
is not related in a simple way to the local parameters of the 
gaseous disk. Our simulations show that the simple reconstruc- 
tion of the scale height of the gaseous disk is possible only from 
distribution of giant planets at greater distances than ~ 8 AU, 
as planets at such large orbits did not significantly move away 
from their original places. 

If we include in our models that the planet can still grow in 
mass after gap opening, it is even more difficult, if at all possi- 
ble, to derive conclusions about the gaseous disk (its tempera- 
ture, scale height, etc.) based on the distribution of giant plan- 
ets. This is because the rate of this additional accretion onto the 
planet depends mainly on the accretion rate in the disk, while 
the time frame in which this accretion can occur is determined 
by the surface density of the planetesimal swarm. The effect 
of the additional growth of the planet can even lead to rever- 
sion of the correlation between the masses of the planets and 
semimajor axes of their orbits. 

We also checked how the mass of the central star influ- 
ences the distribution of planets. The main result is that in disks 
around less massive stars, giant planets at the given location 
tend to be less massive. At the same time, the giant planets 
with a given mass tend to form closer to the less massive stars. 

Our models ca n be seen as an exte nsio n of previous calcula - 
tions performed by Ida & Lin (2004a) and Alibe rt et alj J2005). 
In comparison, the surface density of the planetesimal swarms 
in our models are computed self-consistently. This additional 
factor makes it much more difficult to predict the structure of 



the gaseous disk from the distribution of giant planets much 
more difficult. Nevertheless it does seem necessary. Whenever 
we performed the similar calculations as described above, but 
with the assumption that the surface density of the planetesimal 
swarm is always a constant fraction of the surface density of 
gas, we were not able to perceive any giant planets. This seems 
contradict the results of those two papers. A possible reason 
for this different result may be the difference in the underlying 
models of the gaseous disk. However, we think that the main 
factors are the different times for the dispersion of the gaseous 
disk (3 Myr vs. 10 Myr) and the different sizes of solids at 
time (small grains vs. planetesimal sizes). Our models show 
that the time needed to reach planetesimal sizes can be as long 
1 Myr and is not negligible. 

Our results also seem to contradic observations that show 
that more massive planets tend to be located farther away from 
the host star. This correlation i s reproduced better by the cal- 
culations of llda&Lirl |2Q04a), who conclude that maximal 
mass of the planets after including type II migration does not 
strongly depend on the size of the o rbit. On the other hand, in 
the results from Al ibert et alJ J2005t) the correlation is identical 
to ours, namely the more massive planets tend to be closer to 
the star . This difference can have two sources. First. Ilda & Linl 
(2004a) neglect the accretion of gas onto the planet, after it 
opens a gap in the disk. The second reason, which we think 
is more important, can again be the difference in the structure 
of the gaseous disk. In our calculations ratio (H/r) is mainly a 
decreasing function of the radius in the inner parts of the disk, 
while it is an increasing function in the disk model of Ida & Linl 
(2004a). Because this ratio determines the mass of the planet 
which in turn is able to open a gap in the disk, then it can play a 
big role in determining the properties of the set of giant planets. 
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